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It was recently proposed that the electron- frame dissipation measure, the energy transfer from the electro- 
magnetic field to plasmas in the electron's rest frame, identifies the dissipation region of collisionless magnetic 
reconnection [Zenitani et at. Phys. Rev. Lett. 106, 195003 (2011)]. The measure is further applied to the 
electron-scale structures of antiparallel reconnection, by using two-dimensional particle-in-cell simulations. 
The size of the central dissipation region is controlled by the electron-ion mass ratio, suggesting that electron 
physics is essential. A narrow electron jet extends along the outfiow direction until it reaches an electron 
shock. The jet region appears to be anti-dissipative. At the shock, electron heating is relevant to a magnetic 
cavity signature. The results are summarized to a unified picture of the single dissipation region in a Hall 
magnetic geometry. 



I. INTRODUCTION 

Collisionless magnetic reconnection drives explosive 
events in space plasma environments. On the large scale, 
reconnection is an MHD-scale process but is facilitated 
by a compact "diffusion region" surrounding the recon- 
nection point, where kinetic physics plays a role. In 
the Hall reconnection model, ^ ^ it is thought that the 
electron diffusion region (EDR), where the electrons de- 
couple from field lines, is embedded in the ion's diffu- 
sion region^. The EDR is related to a localized dissi- 
pation region that allows fast reconnection, and it also 
adjusts the overall system evolution to the outer ion 
physics. Owing to its importance in reconnection, the 
EDR is one of the most important targets of NASA's 
upcoming Magnetospheric MultiScale (MMS) mission 
(http://iiims.space.swri.edu/), which will be the first 
to probe the electron-scale physics. 

The Hall reconnection model has been recently chal- 
lenged by kinetic particle-in-cell (PIC) simulations that 
sufficiently resolve electron-scale structureJ^M^. It has 
been argued that the EDR is no longer localized but it 
unexpectedly stretches in the outflow directionJ^. Fur- 
thermore, it appears to embody a two-scale structure of 
inner and outer region Minor differences aside, the 
inner region is a compact region containing the reconnec- 
tion site. It features a strong out-of-plane electron cur- 
rent and a dissipative electric fleld E' = E-\-Ve x S 7^ to 
transport the magnetic flelds. The outer region extends 
in the outflow direction, accompanied by a fast electron 
jet, often denoted as the "super- Alfvenic" jet. Since the 
jet outruns the moving magnetic flelcP, the out-of-plane 
component of E' has the opposite polarity from that in 
the inner region. Satellite observations corroborate these 
structures in near-Earth reconnection site^^^. 
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The role of the EDR is a subject of recent debate. In 
principle, many works have agreed that the reconnection 
rate is controlled by the inner region or the inn er reg ion 
adjusts its size to the global reconnection rat d^ l ^ l -'^-^ -^t^. 
Its long-term, time-dependent behavior remains unclear. 
Daughton et al. ^ argued that the (inner) EDR stretches 
in the outflow direction. They further proposed an idea 
that the secondary island formation maintains fast recon- 
nection by cutting and shortening the elongated EDR. In 
contrast. Shay et al. ^ demonstrated that the inner region 
remains at the flnite length and that the reconnection is 
fast, quasi- steady, and laminar. Regarding the outer re- 
gio n, it does not appear to constrain the reconnection 
rat^^l^^KISI^ because reconnection remains fast while the 
electron jet extends a large distance. Hesse et al.l^ found 
that the electron flow essentially consists of the E x B 
convection and the diamagnetic current for the fleld re- 
versal. Thus, even though it is frequently mentioned by 
the outer EDR^, it is not clear whether the outer region 
is dissipative. Furthermore, it is not well understood how 
the electron jet is terminated. 

Recently, Zenitani et al.^ proposed a new measure to 
identify the dissipation region of reconnection. Consider- 
ing the energy dissipation in the rest frame of electron's 
bulk motion, they introduced the electron-frame dissipa- 
tion measure^ 

De = 7e [j -{E^VeXB)- p,{Ve ' E)] , (1) 

where 7e = [1 — ('Ug/c)^]"^^^ is the Lorentz factor for 
the electron velocity and pc is the charge density. This 
measure is a Lorentz- invariant, and it is related to the 
nonideal energy conversion, which is essential for recon- 
nection problems. Their PIC simulations demonstrated 
that De identifles a single dissipation region surrounding 
the reconnection site. This is substantially different from 
the previous two-scale picture. 

The purpose of this paper is to organize our knowledge 
on the dissipation region in collisionless magnetic recon- 
nection in line of recent theoretical progres^^^. Us- 
ing high-resolution two-dimensional PIC simulations, we 
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perform a detailed investigation of internal structures in 
basic antiparallel configurations. In particular, the previ- 
ous two-scale EDR are reexamined by using the electron- 
frame dissipation measure D^. The structures are better 
understood in a Hall magnetic geometry. We further ex- 
plore the termination region of the electron jet. 

This paper is organized as follows. We briefly review 
the dissipation measure in Sec. [llj We describe our 
numerical setup in Sec. |III[ We present the simulation 
results in Sec. IV We carefully investigate fine struc- 
tures in the following subsections. Section [V| contains 
discussion and summary. 



II. THE ELECTRON-FRAME DISSIPATION MEASURE 

For better understanding, we review a basic concept 
in Ref. 16 within the nonrelativistic physics. The essence 
is to evaluate the Ohmic dissipation in the rest frame 
of electron's bulk motion. Let the prime sign C) denote 
quantities in the electron frame. A simple algebra yields 



erijV. 



PcVe- 



e{niVi - UeVe) - e{ni - ne)Ve 



(2) 



The last term comes from the convection current (pcVe)^ 
arising from the motion of the charged frame. The elec- 
tric field in the electron frame is given by the nonideal 
electric field, 

= E^VeX B. (3) 
The energy dissipation in the electron frame yields 

De = f -E' =j-{E^VeXB)- p,Ve ' E. (4) 

This is an invariant scalar with respect to the Galilean 
transformation and is equivalent to Eq. [l]in the limit of 
7e ^ 1. We employ this nonrelativistic formula through- 
out this paper. We often discuss the composition of 
in the following form. 



De = j:rE'^ + jyE'y + j,E'^ - p^Ve ■ E. 



(5) 



The last term deals with the effective energy conversion 
by the convection current. Hereafter we refer to it by the 
"charge term." 

In a neutral plasma, one can obtain = j ■ E' by 
dropping the charge term. This reduced form is equiva- 
lent to the Joule dissipation in the plasma rest frame in 
Birn & Hesse^f^, who discussed the energy transfer in 
MHD and kinetic reconnection systems. 



III. NUMERICAL METHOD 



We use a partially-implicit PIC code^^E^. The 
length, time, and velocity are normalized by the 
ion inert ial length di = c/upi^ the ion cyclotron 
time = mic/{eBo), and the ion Alfven speed 



CAi = 5o/(47rmino)"^/^, respectively. Here, ujpi = 
(47re^no/m^)^/^ is the ion plasma frequency, no is the 
reference density, and is the reference magnetic 
field. We employ a Harris-like configuration, B{z) = 
Botdinh.{z/L)x and n{z) = no[0.2 + cosh~^(z/L)], where 
L = 0.5di is the half thickness of the current sheet. The 
electron- ion temperature ratio is T^/Ti = 0.2. Periodic 
{x) and reflecting wall {z) boundaries are used. 

In the main run (run lA), the mass ratio is mi/rrie — 
100. The ratio of the electron plasma frequency to the 
electron cyclotron frequency is uOpe/^ce = 4. In this case, 
the speed of light is {c/caz) = {mi/meY^'^{ujpe/Vtce) = 
40. The domain of [0, 76.8] x [-19.2, 19.2] is resolved by 
2400 X 1600 cells, and 2.2x10^ particles are used. We 
carry out several other runs, changing the two key pa- 
rameters: the mass ratio m^/mg and the uOpe/Vtce Pa- 
rameter. All parameters are presented in Table, [ij In 
all cases, reconnection is triggered by a small flux per- 
turbation. 8Ay = -2L5iexp[-(x^ + z^)/(2L)^], where 
Bi = 0.1^0 is the typical amplitude of the perturbed 
fields. The initial current is set up accordingly. 



TABLE I. Simulation parameters 



Run 


rrii/me 


UJpe 1 ^ce 


Domain size 


Grid cells 


Particles 


lA 


100 


4 


76.8 x 38.4 


2400 X 1600 


2.2 X 10^ 


IB 


100 


2 


76.8 X 38.4 


2400 X 1600 


2.2 X 10^ 


2A 


25 


4 


102.4 X 51.2 


1600 X 1600 


1.4 X 10^ 


2B 


25 


2 


102.4 X 51.2 


1600 X 1600 


2.6 X 10^ 



IV. RESULTS 



Overview 



We identify the dominant reconnection site by find- 
ing the minimum in x of the magnetic flux <l>(x,t) = 
\l\Bx\dz. Then we discuss the system evolution by us- 
ing the normalized flux transfer rate or the reconnection 
rate 7^, 



7^ = - 



(6) 



where the subscript in denotes quantities at 3di upstream 
of the dominant reconnection site. 

Figure [l] shows the time evolution of 71. The rate 
quickly increases, modestly overshoots at t~15, and then 
approaches a quasisteady value of 0.12-0.13. Such qua- 
sisteady evolution corroborates previous investigations.!^ 
The rates are insensitive to the mass ratio in this 
normalization.!^ The parameter cj^e/^ce does not make 
a significant difference also. We confirmed that the re- 
connection rates are independent of domain sizes in the 
timescale of our interest, by carrying out supplemental 
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FIG. 1. The normalized reconnection rates IZ (Eq. [g]) as a 
function of a time (Qcit). 



small runs. Meanwhile, as will be discussed later, the pe- 
riodic domain effects start to modulate the outflow struc- 
tures after t~40 in runs lA and IB. 




FIG. 2. (Color online) Magnetic field line structure in run 
lA, averaged over t = 35-36. Rear panel: the out-of-plane 
magnetic field By . Front panel: the electron-frame dissipation 
measure De 



Figure |2] presents a snapshot of run lA. Physical quan- 
tities are averaged over t = 35-36. At this stage, the re- 
connection is steadily going on as can be seen in Figure [l] 
and the characteristic structures are well developed. This 
is the stage of our primary interest. The gray lines show 
magnetic field lines in 3D. Magnetic field lines are flipped, 
because they tend to follow the electron's flow in the —y 
direction. This is a well-known signature of Hall recon- 
nection. For example, Huba & Rumanov*^ showed a sim- 
ilar field-line structure by using Hall MHD simulations. 
For example, Yamada et al.^^ found a similar structure 
in an experiment device. Our fully-kinetic simulation is 
consistent with these results. The rear panel in Figure 
[2] shows the out-of-plane magnetic field By . The color 
indicates the polarity: By > in red and By < in blue. 
The projection of flipped fi eld lin es leads to a charac- 
teristic "quadrupole" patterrP^^ in By. In addition, as 
indicated in Figure [2j there are two "magnetic cavities" 



in the outflow region, where By becomes weaker. This 
signature is often found in PIC simulations of Hall recon- 
nection, but this has never been discussed before. The 
front panel shows the electron-frame dissipation measure 
De (Eq. |4|. A central red region of I^g > is the dissi- 
pation regioiP^. We will investigate this region in detail 
later in this paper. 

Figure [3] show various averaged quantities over t = 35- 
36 in run lA. Figure [sj^a) shows the electron outflow ve- 
locity Vex- One can see bi-directional jets from the recon- 
nection site. They travel much faster than the upstream 
Alfven speed (0.2)~^/^c^i ~ 2.24, which approximates 
the outflow speed in an MHD scale. The electron jets 
are remarkably narrow. Their widths are on an order of 
the local electron inertial length. Outside the separatri- 
ces, there are weak reverse flows toward the reconnection 
site. In Figure |3jb) we show the electron current density 
|je|. It is distinctly strong in a narrow region near the 
neutral plane {z = 0). We note that it is remarkably 
flat in the x-direction. Thus, \je\ appears to be a good 
marker of the interesting region. Hereafter we call the 
region an "electron current layer." 

Figure |3|c) shows the vertical electric field Ez. Near 
the reconnection site, one can see a bipolar signature 
across the electron current layer: is negative for z > 
and positive for z < 0. This is called the Hall or polar- 
ization electric fielcP'^^M^. This is supported by charge 
separation between a broader distribution of meandering 
ions and a narrow distribution of meandering electrons. 
The bipolar peaks \Ez\ are typically 8 times stronger 
than the reconnection electric field Ey. In the inflow re- 
gion, outside the electron current layer, magnetized elec- 
trons tend to travel in the -y-direction^^. In the outflow 
region, the Ez region consists of two parts: the large- 
scale X-shaped structure along the separatrices and the 
top/bottom boundaries of the electron current layeif^. 
We think the latter is related to the electron meandering 
motion, similarly as the central region. In addition, as 
mentioned by a previous work^^, one can recognize an 
inverted signature of Ez in a very vicinity of the neutral 
plane {\z\ < 0.1) inside the electron current layer. For 
example, Ez is positive for z > and negative for z < 
around x ^ 29 and x ~ 47. These "inversion electric 
field" are supported by charge separation inside the 
electron meandering orbit. 

Shown in Figure [sj^d) is the out-of-plane component 
of the electron Ohm's law, Ey = [E -\- Ve x B]y^ in unit 
of caiBq. One can recognize two distinct regions along 
the electron jets. One is the compact region surrounding 
the reconnection site, where Ey is positive. The other re- 
gions extend in the outflow direction and Ey is negative, 
because the fast electron jets outrun the field convec- 
tion. Those two are often referred as the inner and outer 
EDRs in recent literaturP^l. Note that E'y ^ just tells 
us that the ideal assumption breaks down. It does not 
always indicate magnetic dissipation. 

Figure [sj^e) presents the electron-frame dissipation 
measure (Eq. [4| in units of caiBqJq^ where Jo is the 
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FIG. 3. (Color online) Averaged profile of run lA over t = 35-36. (a) The electron outflow velocity Vex^ (b) the electron current 
density \je\ in units of Jo, (c) the electric field Ez in units of CAiBo, (d) the out-of-plane component of the electron nonideal 
condition {E-\-Ve x B)y in units of CAiBo, (e) the electron-frame dissipation measure De in units of CAiBoJo, and (f) the charge 
term {—pcVe • E) in De. The contour lines are in-plane magnetic field lines. The dash line indicates the field reversal, Bx — 0. 



initial Harris current density. Note that the color scale 
is different from a reduced color scale in Fig. [2] One can 
see that Dq gives a different picture from Figure [sj^d) . 
The dissipation region with > {) \s located near the 
reconnection point [the red region near 35 < x < 41 in 
Fig. |3je)]. On the other hand, is weakly negative 
in the outflow regions, as can be seen in light green in 
the Figure. There, its amplitude \De\ is substantially 
smaller than near the reconnection point. Figure [sj^f) 
presents the last charge term in Eq. |5]in the same unit 
as Fig. [sj^d). This term appears only in the close vicinity 
of the above central region of I^g > 0? and it is an order 
of magnitude smaller than D^. 



B. Dissipation region 

Here, we focus on the dissipation region. The panels 
in Figure [4] present the electron nonidealness E'y^ the dis- 
sipation measure D^^ and the dissipation measure in 
run 2A at lower mass ratio of mi/rrie = 25. We find that 
the lengths of the electron current layer and the dissi- 



pation region are nearly stationary in the later phase of 
simulations. The panels in Figure |4] are taken from the 
stationary phase. 

Comparing Figs. [4]^ a) and (b), one can see that the 
dissipation region is longer in t he o utflow direction than 
the inner region defined by E'^lIU^. The comparison of 
Figs. [4jb) and (c) suggests that the size of the dissipa- 
tion region is controlled by the electron mass. In our 
ion-based units, the electron inertial length scales with 
(me/mi)^/^ and the electron bounce width scales with 
[me/miY^^ (Refs. TOlEOj). We find that the length and 
the width of the dissipation region in run lA is smaller 
than those in run 2 A by a factor of 1.5-1.8. This is in- 
termediate between (100/25)^/"^ and (100/25)^/^ The 
scaling to the mass ratio needs further investigation, but 
electron physics appears to be essential for the dissipa- 
tion region. 

Next, we discuss the composition of the dissipation 
measure D^. Shown in Figure [5] are contributions from 
the terms in Eq. [s] and the charge separation (n^ — 
^e)/(^i + ^e) along the inflow and outflow lines. The 
shaded regions indicate the dissipation region by the sign 
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(a) [E+v,xB] 




FIG. 4. (Color online) (a) the out-of-plane component of the 
electron nonideal condition (E -\- Ve x B)y around the recon- 
nection site, (b) the electron-frame dissipation measure De, 
and (c) De in run 2A. 



{De > 0; orange) and the e- folding length of (the in- 
ner red region), respectively. The blue dotted lines indi- 
cate the electron ideal condition, Ey = [E-\-VeX B]y =0. 
The nonideal electric field Ey is positive between the two 
dotted lines. 

In the inflow direction in run lA [Fig. (Sja)], De is well 
localized near the neutral plane z = 0. Its main contrib- 
utor is the y-teim (jyEy). The charge term {—pcVe-E) is 
responsible for —4% of D^- This term has a double peak 
structure and comes from a double-peak charge distri- 
bution (the red dashed line), due to the electron bounce 
motion. The charge separation is up to —4% in this case. 
The z-term {jzE'^) has also double peaks, because both 
jz and E'^ have bipolar signatures. However, since jz is 
small, the z-term is a very minor contributor to D^. The 
x-term is negligible also. 

Figure [5|b) shows the composition in run IB at lower 
ratio of uj^^l^ce — 2. In this case, the charge term is 
responsible for -15^20% of D^^ because the charge sepa- 
ration is significant, ~ 10%. A lower uOpe/Q^ce allows more 
significant charge separation. We confirm this trend at 
the mass ratio mi/rrie — 25, too. In the outflow direction 
[Figure [sj^c)], one can see that the x-term {^jx^'x) 
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FIG. 5. (Color online) Composition of the dissipation measure 
De (a) along the inflow line [x — 37.9) in run lA, (b) along 
the inflow line (x — 37.3) in run IB, and (c) along the outflow 
line (z = 0) in run lA, 



other contributor to D^- Both D^ and ix^'x noisy in 
the outflow regions, because Ex is noisy. The gray region 
is defined by the spatial location of a maximum electron 
outflow velocity. 

Let us compare the dissipation region with the conven- 
tional inner region^ ^. In the inflow direction [Fig. [sja)], 
the inner region has been defined by the electron non- 
idealness E^ ^ (the blue dotted lines). Both orange 
and red dissipation regions are thinner than the conven- 
tional inner region, bounded by the blue dotted lines [Fig. 
(Sja)]. One case see that the e-folding region (red) is rele- 
vant to the electron bouncing region, featuring the charge 
separation. On the other hand, there have been several 
definitions in the outflow direction: the electron outflow 
velocity (gray region j^^ ^, the sign of the electron nonide- 
alness (the dotted lines)P^^I^, or the out-of-plane current 
profile^^. Note that Klimas et al.l^ employed a different 

formula — [(l/nee)V • ^ + {me/e){ve ■ ^)ve)]y > 0, but 
this is almost equivalent to the popular criteria Ey > 0. 

In the outflow direction, the dissipation region is 
longer than the inner region defined by ^^HU^. One 
can intuitively understand this, considering the magnetic 
geometr}/^^. Since the field lines are flipped (Fig. we 
consider the curved surface as illustrated in Figure [6] 
Hereafter we call it the Hall surface. As the field lines 
rotate their directions in the surface, j (^ V x B) is nor- 
mal to it. The J -aligned (out-of-plane) component of 
is responsible for the nonideal field convection in the Hall 
surface. Meanwhile, at sufficiently high Upe/^ce^ one can 
reduce the dissipation measure to D^ ^ j • E' . This is 
controlled by the j-aligned component of E' . In other 
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FIG. 6. (Color online) The Hall structure of the dissipation 
region. 



words, De autoDiatically takes care of the out-of-plane 
component of E' with respect to the Hall surface. Thus, 
as we move to the outflow direction, we should consider 
E'^ as well as Ey. The positive- I^e dissipation region is 
always longer than the positive- regioiP^, because 
further includes contributions from E'^. According to our 
preliminary results (Klimas et al. in prep), the former is 
commonly twice longer than the latter, even when the 
reconnection region is disturbed by secondary magnetic 
islands. 

On the other hand, the dissipation region is shorter 
than the region defined by the peak location of Ve^^- 
We find that v^x is not very useful to define an important 
region due to the following reasons. First, it is difficult 
to identify the peak location, when gradually changes 
in the outfiow direction. Second, the physical meaning of 
the maximum v^x is unclear in the Hall geometry (Fig. 
[6|. The Hall surface is already rotated around there, and 
so. Vex is highly infiuenced by the out-of-plane diamag- 
netic fiow^ and the angle of the Hall surface. 

Practically, we prefer to define the dissipation region 
by the e-folding length (red) rather than by its sign (or- 
ange). This is because in some cases the slope of 
is very fiat around ^ {) m the infiow directions. In 
Figures |5|a) and (c), the aspect ratios of the dissipation 
region are 0.39 : 4.77 (= 0.08 : 1) by the e-folding length 
and 0.77 : 6.83 (= 0.11 : 1) by the positive dissipation, 
respectively. These values reasonably fit to the typical 
reconnection rate of 0.1. 



C. Outflow region 

Next, we visit the outfiow regions [29 < x < 36 or 
41 < X < 48; Fig. |3ja)]. One can see the narrow fast 
electron jet They are fianked by slow uniform fiow 
regions. The region is denoted as the pedestaF^. In 
the pedestal, the Ey component is negligible, suggesting 
quasi-ideal electron fiows. There is a boundary layer with 
weak E'^ and E'^ between the jet and the pedestal. 

Importantly, as shown in Figure [sj^e), the dissipation 
measure is weakly negative < {) m. the electron jet. 
From the viewpoint of the Hall surface (Fig. [6|, the out- 



of-plane component (~ j-aligned component) of E' is not 
exactly zero in a rotated frame^^. The electron jet still 
outruns the field convection in the outfiow direction in 
the Hall surface. This leads to ^ j ■ E' < in the 
electron jet. On the other hand, is weakly positive 
in the pedestal and it often has positive peaks in the 
boundary layers between the jet and the pedestal. This 
is visible in particular in the downstream (29 < x < 33 
and 45 < X < 48) and is more evident at later stages. 

The fact De < tells us that the electron jet region 
is not dissipative. To understand this, it is useful to 
discuss the MHD energy balance, as given in previous 
liter at ur^I^KI^. We consider a mass- averaged MHD ve- 
locity in an electron-ion plasma. 



'^mhd = 



miTliVi + meUeVe 



(7) 



The energy transfer from the fields to plasmas in the 
MHD frame I^mhd is similarly given by 

^mhd =j-{E^ Vmhd X B)- Pc{Vmhd * E) . (8) 

A small amount of algebra leads to 



{rrii + me)ne ^ 

^mhd = \ J^e- 



(9) 



Thus De < indicates l^mhd < i-e., the plasma's en- 
ergy is converted to the field energy even in the MHD 
frame. This contradicts with the concept of magnetic 
dissipation or diffusion, which consumes the field energy 
in the MHD frame. No dissipation process takes place 
here. Instead, this region is better described as "anti- 
dissipative." 

Figure [TJ^a) shows a flip angle of Hall magnetic field 
lines. We calculate it by an arctangent of By over Bx, 
which is a good indicator outside the z = plane. The 
sign of the angle is neglected to better see the structure. 
As can be seen in Figure [2j the magnetic field lines are 
extremely fiipped. The angle is almost 90° in the elec- 
tron jet region. Such a strong fiipping is found in recent 
studies (e.g., the hodogram analysis by Drake et al.^^. 
The angle changes from zero to 90° from the upstream 
region to the electron jet region. Further observation tells 
us that the angle has a two-scale structure in the same 
x-location: a moderately fiipped pedestal and a highly 
fiipped electron jet. This is more evident at a later time 
of t = 40-41 [Fig. [7|b)] and the ID cuts along the ver- 
tical lines [Fig. [tJ^c)]. In the pedestal region, the angle 
is 45 — 50°. Note that the Hall magnetic field \By\ is 
usually comparable to the antiparallel field l^^^l in many 
simulations. 

We think the Z^g-pattern and the field line geometry is 
relevant. In the electron jet, the field lines are strongly 
fiipped, and < implies that the electromagnetic 
fields gain energy in the electron frame. Outside the jet, 
the magnetic field lines are mildly fiipped and > 
indicates that plasmas gain energy from the fields. The 
electron-field energy transfer • E shows similar signa- 
tures as De. We speculate that the bulk kinetic energy 
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(a) Flip angle 




angle 



FIG. 7. (Color online) The flip angle [degeree] of the Hall 
magnetic field lines (a) at t = 35-36 and (b) at t = 40-41. (c) 
ID cuts along the two lines. 




FIG. 8. (Color online) (a) The stack plot of the electron out- 
flow speed Vex (thin lines) and the reconnected magnetic field 
Bz (thick lines), (b-e) Averaged properties near the electron 
jet front at t = 35-36: (b) Bx, (c) the out-of-plane magnetic 
field By, (d) the electron temperature ^{pxxe +Pyye +Pzze) 
in units of me(?Ai^ C^) ^he electron energy transfer je ■ E 
in units of CAiBoJo. 



of the fast electron jet is transferred from the jet to the 
pedestal via the magnetic field lines that thread two re- 
gions. This process needs to be further investigated. 



D. Electron shock 

We find a new shock-like structure at the fast electron 
jet front. Shown in Figure [sj^a) are the time evolution of 
the electron outflow velocity (vex) and the reconnected 
magnetic field (Bz). At t = 36, one can see that Bz be- 
comes twice stronger (from ~0.08 to '^0.16) across the 
transition region, x ~ 49. On the other hand, the fast 
electron jet suddenly slows down there, and then the 
electrons are magnetized further downstream. We also 
find that the electron temperature and pressure increase 
across the transition region. These signatures indicate 
that it is a shock between the upstream fast electron jet 
and the downstream magnetized electron flow. The dis- 
continuity separates an unmagnetized upstream flow and 
the magnetized downstream flow. Ions are insensitive to 
the shock. They are not magnetized on both sides, and 
the length of the transition region is only one (i^, which 
is smaller than the ion inertial length based on the lo- 
cal density. In order to distinguish this from full MHD 
shocks, we hereafter refer to it as an "electron shock." 



As time goes on, the electron shock moves outward. 
It reaches at x ~ 53 at t = 44. After that, the shock 
propagates backward. This is the influence of our peri- 
odic boundary condition. After a sufficiently long time, 
the back pressure from the downstream region is strong 
enough to push the shock front backward, like a reverse 
shock reflected by a wall. The shock becomes stronger 
than the earlier phase, as seen in B^ in Figure [sja). 

Shown in Figures[8]^b)-[8je) are physical quantities near 
the electron shock region. Interestingly, B^ changes its 
sign in the shock transition region of 48 < x < 49: 
< on the upper side and > on the lower side 
[Fig. [sj^b)], due to a reverse electron current jey < 0. 
Unlike the dissipation region, electrons travel in the -\-y- 
direction there, because the fast-traveling electrons are 
just trapped by B^. In the upstream region of x < 48, By 
sharply changes its polarity near the electron jet, while it 
is rather uniform outside the jet [Fig. [sj^c)]. In the shock- 
downstream of X > 48, By suddenly becomes weaker 
than in the upstream. The vertical magnetic pressure 
^{B'^ -|- By) exhibit similar profiles. This is identical to 
the magnetic cavity (Figure |2|. Note that the vertical 
component {Bz) is compressed across the shock. 

We calculate the electron temperature [Fig. [sj^d)] by 
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tracing the diagonal components of the pressure tensor, 



— Q {Pxxe ~^ Pyye ~^Pzze)- 



(10) 



Each component basicahy shows a similar profile. One 
can see that the temperature quickly increases across the 
shock, due to the fast jet speed. The energy transfer 
from the electromagnetic fields to electrons, • is 
presented in Figure [sj^e) . It is negative je - E < around 
the transition region. The same signature of je • E < 
was reported by a previous work (the "electron dynamo 
region" in Ref. l3T]) . 

These results tell us that the bulk kinetic energy of 
the electron jet is transferred to the electron heat [Fig. 
jsjd)] and the magnetic energy [Fig. [sj^e)] in the shock 
downstream. Here, the magnetic energy is stored in the 
compressed reconnected field 5^, because both and 
By are very weak near the neutral line z = 0. The strong 
electron heating also explains the magnetic cavity struc- 
ture. Although the electron pressure substantially in- 
creases across the shock, the surrounding tangential fields 
{Bx^By) are not strong enough to confine electrons in the 
z-direction in the downstream. The high-pressure elec- 
trons expand the structure outwards along Bz in the ±z 
directions. The displacement of the tangential fields lead 
to the magnetic cavity. 



V. DISCUSSION 

Starting with a work by Daughton et al. ^, the struc- 
tures of the electron diffusion region has been actively 
discussed for the past five years. Earlier investigations 
focused on the violation of the ideal frozen-in condition, 
E' ^ {). This lead to the popular two-scale pictur^^, 
in which the role of the EDRs was not clearly under- 
stood. In this work, we have explored and reorganized 
our understanding of fine reconnection structures in line 
of Ref. l16.. We have confirmed that the electron-frame 
dissipation measure characterizes the critical region at 
higher mass ratio of m^/mg = 100. Since the size of the 
dissipation region becomes smaller at higher mass ratio, 
it seems that the electron physics is responsible for the 
dissipation region. The relevance to the mass ratio fur- 
ther needs to be investigated in order to extrapolate our 



results to a realistic ratio of 1836. In Sec. IV we found 



that no dissipation takes place in the electron jet region, 
^mhd < 0. Consequently, we don't think it is appropriate 
to call it the outer electron "diffusion" region or the outer 
electron "dissipation" region. As reconnection proceeds, 
the nondissipative electron jets are elongated, while the 
central dissipation region remains compact. This is con- 
sistent wi th t he recent consensus that reconnection re- 
mams 

From the observational viewpoint, the contribution 
from the charge term in Eq. [5] is of strong interest, be- 
cause the charge density pc will be difficult to probe. In 
run lA, the charge term is an order-of-magnitude smaller 



than Dq^ and it is only localized near the dissipation 
region [Fig. (sjf)]. In run IB with ujpe/^ce = 2, the 
charge separation appears nonnegligible. Let us estimate 
the charge separation around the dissipation region. We 
assume that the local plasma density 2n = + ng is 
roughly uniform. Electrons are magnetized to the field 
lines outside the electron current layer, and the electron 
out-of-plane speed is on an order of the electron Alfven 
speed CAe- The Hall electric field Ez is approximated by 



Ez 



Vey\BQ (CAe 



Gauss's law in the z direction tells us the charge density 
in the center. 



e(ni - He) 



(12) 



where 6 is the typical thickness of the dissipation region. 
Using the result in Sec. |IVB[ we assume 



(-)(-)^' 



(13) 



where the index s is between 1/4 < s < 1/2. Substituting 
this to Eq. [121 we obtain 



rij 



\meJ 



1 

Siren \uj, 



pe 



We see that the charge separation is controlled by the 
Upe/^ce parameter and that it is much less sensitive to 
the mass ratio. 

In the magnetotail, the magnetic field is Bq ^ 20 nT 
in the lobe (upstream region). The corresponding elec- 
tron gyro frequency is flee ^ 560 Hz. The plasma sheet 
density no ^ 1 cm~^ gives Upe ^ 9 kHz. Thus we 
obtain a typical value of uOpe/Vtce ^ 16. In addition, 
since the upstream plasma occupies the reconnection re- 
gion, it is useful to compare the upstream conditions. 
While our numerical model assumes that the background 
density is 0.2no, the lobe density is substantially lower, 
O(O.l) - (9(0.01) cm-^. In the case of 0.02 cm"^, since 
an upstream plasma frequency is modified by a factor of 
aAO, our model with cj^g/^ce ^ 5 will be relevant. These 
estimates suggest that the charge separation effects will 
be marginally negligible in the magnetotail. This needs 
to be verified by further numerical investigations with 
realistic lobe densities^^ 

We have found that the fast electron jet terminates at 
the shock, where the electrons become magnetized. The 
magnetic cavity structure can be understood as its con- 
sequence. We expect that shock-like signatures are more 
prominent at higher mass ratio, because the electron jet 
speed ^CAe is much faster than the typical outflow s peed 
^CAi- In a larger system or an open-boundary system^^, 
we expect that the electron shock travels further down- 
stream, unless (1) an obstacle in the downstream such 
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I (G) Electron shock 
J & magnetic cavity 



(D) Dissipation region 
(C) Electron current layer (e) Electron diamagnetic jet 



FIG. 9. (Color online) Our present understanding of Hall reconnection structure: (A) Quadrupole magnetic field By fRefs. l22l 
and 12 31 , (B) Hall current system (Ref. [22}, (C) electron current layer (Refs. [4] and (5}, (D) dissipation region fRef. [TBI Sec. 
IV B), (E) electron diamagnetic jet fRefs. [6l [7[ [T5l Sec. IV C), (F) pedestal (Ref. TT'), and (G) electron shock and magnetic 
cavity (Sec. |IVD|). 



as a secondary island, (2) instabilities of the jet, and (3) 
3D effects interrupt the shock propagation. The electron 
shock may play a role on particle acceleration as well. For 
example, Hoshino et al. discussed a two-step scenario 
of electron acceleration: pre-acceleration near the recon- 
nection site and energetization by VB/curvature drifts 
near the flux the flux pile-up region. Our results suggest 
that a significant shock-heating (and possibly shock-drift 
type acceleration) takes place between the two accelera- 
tion sites. The electron shock deserves further investiga- 
tion in the context of particle acceleration. 

We summarize our understanding in Figure [21 As well 
known, a Hall reconnection features the quadrupole out- 
of-plane magnetic field By , which is generated by the Hall 
current circuit^^. There is a narrow channel of Ey ^ 
near the neutral plane^. We call this channel the elec- 
tron current layer, because the electric current \je\ is dis- 
tinctly stronger than in the other regions. This current 
layer is beyond the scope of the fluid theory. Inside the 
electron current layer, there is a single dissipation re- 
gion surrounding the reconnection point^^. This region 
is thinner and longer than the conventional inner region. 
The fast, narrow electron jeiP^ travels from the dissi- 
pation region. This is a projection of the diamagnetic 
current and a part of the Hall current circuit—. As 
discussed in Section [IV C [ this region is not dissipative. 
There is a pedestal region''^^ outside the electron jet. The 
electron jet is terminated by an electron shock and the 
magnetic cavity develops there. 

The upcoming MMS mission will measure the field and 
plasma properties at high resolutions in near-Earth re- 
connection sites. Due to the limited bandwidth between 
the satellite and the Earth, it is very important to se- 
lect the data in high-priority regions. We have demon- 
strated that the electron-frame dissipation measure is a 
good marker of the most important dissipation region. 
We have further updated our understanding on fine struc- 



tures surrounding the dissipation region. We hope that 
the measure is helpful to maximize the scientific profit of 
the mission. 



ACKNOWLEDGMENTS 

The authors acknowledge Keizo Fujimoto and Nicolas 
Aunai for useful comments. One of the authors (S.Z.) 
acknowledges support from JSPS Postdoctoral Fellow- 
ships for Research Abroad. This work was supported by 
NASA's MMS mission. 

ij. Birn, J. F. Drake, M. A. Shay, B. N. Rogers, R. E. Den- 
ton, M. Hesse, M. Kuznetsova, Z. W. Ma, A. Bhattacharjee, 
A. Otto, and P. L. Pritchett, J. Geophys. Res. 106, 3715, 
doi:10.1029/1999JA900449 (2001). 

^M. A. Shay, J. F. Drake, B. N. Rogers, and R. E. Denton, J. 
Geophys. Res. 106, 3759, doi:10.1029/1999JA001007 (2001). 

^J. F. Drake and M. A. Shay, in "Reconnection of Magnetic Fields: 
Magnetohydrodynamics and Collisionless Theory and Observa- 
tions," edited by J. Birn, E. R. Priest (Cambridge University 
Press, Cambridge, 2007), Sec. 3.1. 

^W. Daughton, J. Scudder, and H. Karimabadi, Phys. Plasmas 
13, 072101 (2006). 

^K. Fujimoto, Phys. Plasmas 13, 072904 (2006). 

^H. Karimabadi, W. Daughton, J. Scudder, Geophys. Res. Lett. 
34, L13104, doi:10.1029/2007GL030306 (2007). 

'^M. A. Shay, J. F. Drake, and M. Swisdak, Phys. Rev. Lett. 99, 
155002 (2007). 

8p. L. Pritchett, J. Geophys. Res. 106, 3783, 

doi:10.1029/1999JA001006 (2001). 
^T. D. Phan, J. F. Drake, M. A. Shay, F. S. Mozer, and J. P. 

Eastwood, Phys. Rev. Lett. 99, 255002 (2007). 
i^T. Nagai, 1. Shinohara, M. Fujimoto, A. Matsuoka, Y. 
Saito, and T. Mukai, J. Geophys. Res. 116, A04222, 
doi:10.1029/2010JA016283 (2011). 
iij. F. Drake, M. A. Shay, and M. Swisdak, Phys. Plasmas 15, 
042306 (2008). 

i^A. Klimas, M. Hesse, and S. Zenitani, Phys. Plasmas 15, 082102 
(2008). 

13 W. Wan and G. Lapenta, Phys. Plasmas 15, 102302 (2008). 



10 



^^A. Klimas, M. Hesse, S. Zenitani, and M. Kuznetsova, Phys. 

Plasmas 17, 112904 (2010). 
i^M. Hesse, S. Zenitani, and A. Klimas, Phys. Plasmas 15, 112102 

(2008). 

^^S. Zenitani, M. Hesse, A. Klimas, and M. Kuznetsova, Phys. Rev. 
Lett. 106, 195003 (2011). 

I'^J. Birn and M. Hesse, Ann. Geophys. 23, 3365 (2005). 

isj. Birn and M. Hesse, Ann. Geophys. 27, 1067 (2009). 

^^M. Hesse, K. Schindler, J. Birn, and M. Kuznetsova, Phys. Plas- 
mas 6, 1781 (1999). 

20j. D. Huba and L. I. Rudakov, Phys. Plasmas 9, 4435 (2002). 

21m. Yamada, Y. Ren, H. Ji, J. Breslau, S. Gerhardt, R. Kulsrud, 
A. Kuritsyn, Phys. Plasmas 13, 052119 (2006). 

-^-^B. U. O. Sonnerup, "Magnetic field reconnection," in Solar Sys- 
tem Plasma Physics, edited by L. J. Lanzerotti, C. F. Kennel, 
and E. N. Parker (North-Holland, New York, 1979), vol. 3, p.45. 

23t. Terasawa, Geophys. Res. Lett. 10, 475, 
doi:10.1029/GL010i006p00475 (1983). 



24m. a. Shay, J. F. Drake, R. E. Denton, and D. Biskamp, J. 
Geophys. Res. 103, 9165, doi: 10.1029/97JA03528 (1998). 

25K. Arzner, and M. Scholer, J. Geophys. Res. 106, 3827, 
doi:10.1029/2000JA000179 (2001). 

26p. L. Pritchett, J. Geophys. Res. 106, 25961, 
doi:10.1029/2001JA000016 (2001). 

2^M. Hoshino, T. Mukai, T. Terasawa, and I. Shinohara, J. Geo- 
phys. Res. 106, 25979, doi:10.1029/2001JA900052 (2001). 

28a. Ishizawa, R. Horiuchi, H. Ohtani, Phys. Plasmas 11, 3579 
(2004). 

29L-J. Chen, W. S. Daughton, B. Lefebvre, and R. B. Torbert, 

Phys. Plasmas 18, 012904 (2011). 
30d. Biskamp and K. Schindler, Plasma Phys. 13, 1013 (1971). 
^^M. L Sitnov, M. Swisdak, and A. V. Divin, Phys. Plasmas 114, 

A04202, doi:10.1029/2008JA013980 (2009). 
^^N. Bessho, A. Bhattacharjee, Plasma and Fusion Res. 5, S2017 

(2010). 

33p. Wu, M. A. Shay, T. D. Phan, M. Oieroset, M. Oka, Phys. 
Plasmas 18, 111204 (2011). 



